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Abstract 



^ ' Intermediate mass fragments (IMF) from the interaction of ^''Al, ^^Co and "'^^''Au with 200 MeV 

^ . protons were measured in an angular range from 20 degree to 120 degree in the laboratory system. 

\ The fragments, ranging from isotopes of helium up to isotopes of carbon, were isotopically resolved. 



Double differential cross sections, energy differential cross sections and total cross sections were 
extracted. 

PACS numbers: 
Keywords: 



*'h . machn er @f z-j uelich . de| 

deceased 



1 



I. INTRODUCTION 



Studies of spallation processes, both experimental and theoretical, are numerous. One 
reason for this may be the importance of knowledge of cross sections and reaction mechanisms 
for our understanding of cosmic rays , y, y, y, and the production of cosmogenic 
radionuclides 1^, and the process of neutron production in spallation sources. Recent 
reviews of the process can be found in Refs . Most of the experimental data exist in 

the range above 1 GeV which is important for spallation neutron source construction and 



the understanding of very high energy cosmic rays. However, the energy of the maximum 
abundance of protons in cosmic rays is around 200 MeV 12, ll^ ■ We measured intermediate 
mass fragments at a proton beam energy of 200 MeV incident on three targets spanning 
the periodic table, namely, ^''Al, ^^Co and ^^^Au. These data complement previous cross 
sections for proton, deuteron and tritium emission on ^''Al and ^^''Au jl^ . Q|. The cross 
sections given there for ^He and a-particles are too small when compared to systematics 
[ifil Q| . They were measured with a set up different to that used for the hydrogen isotopes 
and might be low by a factor of 4. We will come back to this point. They also complement 
data for a silver target taken at proton energies close by 
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HQ. 



II. EXPERIMENTS 



The experiment was performed at the separated-sector cyclotron facility of iThemba labs. 
A detailed description of the layout of the facility and equipment is given in Ref. j^] and 
references therein. The beam of 200 MeV was focused to a spot size of less than 2x2 mm at 
the target center of a 1.5 m diameter scattering chamber. Great care was taken to minimize 
the halo of the incident proton beam by focusing the beam through a 3 mm diameter hole in a 
ruby scintillator target. The targets were self supporting foils with thicknesses of 2.9 mg/ cm^, 
1.0 mg/cm^ and 4.0 mg/cm^ for ^^Al, ^^Co and ^^^Au, respectively. The target materials had 
purity of 99.9%. A possible (invisible) oxidation of the surface in case of the aluminum target 
leads to a negligible amount of oxygen. Fragments were measured with a telescope consisting 
of an active collimator followed by three silicon detectors with thicknesses of 50 /im, 150 fim 
and 1 mm. The solid angle of the telescope was 2.2 msr. Another 1 mm thick detector vetoed 
penetrating hydrogen and helium isotopes. The detectors were calibrated with radioactive 
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sources and a precision pulse generator. In order to reduce electronic noise they were cooled 
to a few degrees with chilled water. Detection angles were from 20°to 120°. The opening 
angle of the collimator resulted in an angle uncertainty of ±2.2°. The incident proton flux 
was measured by a beam dump Faraday cup. 

The AE — E method was used for particle identification. A linearized particle identifi- 
cation quantity, PI, was obtained from the energy-range relation, given by 



PI 



{E^ + E^f-El M (1) 



if the particle is stopped in the second detector. Ei denotes the energy deposited in the i-th 
detector and d is the detector thickness. If the particle is stopped in the third detector one 
has the relation 



PI 



{E^ + E2 + E:,f -El /{di + d2). (2) 



Furthermore 



mZ^ oc PI. (3) 



with m the fragment mass and Z the charge number. For the exponent a value b = 1.73 
was used. As an example a mass distribution, obtained by dividing the ranges in the PI- 
spectrum by Z"^, is shown in Fig. ^for the case of cobalt. Hydrogen isotopes fulfilling the 
energy conditions ((H) or Q have the largest yield, but are not considered here. Good isotope 
separation is visible up to boron. In the case of the gold target even carbon fragments could 
be resolved. 

The counting rate was then converted to cross sections. The following systematic errors 
contribute to the total uncertainty. The target thicknesses are known with typically 10% 
uncertainty. The incident flux was measured with 2% uncertainty, while the solid angle, 
electronic dead time correction and energy calibration were estimated to contribute in total 
to less than 2%. The emission angles are uncertain to ±2.2°. The error bars in the figures 
show only the statistical uncertainty. 

In Figs 121- El double differential cross sections are shown for IMFs ranging from '^He to 
^^B. The statistics get poorer with increasing mass number. 

The IMF spectra in the case of the gold target show the effect of the Coulomb barrier: a 
maximum which is in most cases close to 10 MeVxZ with Z the fragment charge number. 
In the case of the cobalt target this is just at the detection threshold which is given by the 
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FIG. 1: Mass distributions from the interaction of 200 MeV protons with ^^Co measured at 20°. 



thickness of the first AE detector. For aluminum the Coulomb barrier is below our detection 
range. In case of the gold target a second component shows up. This is emission below the 
Coulomb barrier of a gold-like system. It is obviously emission from a system with a much 
smaller Coulomb barrier. Unfortunately, the first AE detector is too thick to study such a 
component in the case of the other targets. Such a component can be explained as emission 
from fission fragments which in the case of lighter target nuclei are not as frequent as in the 
case of gold and was also seen in the emission of low energy protons following p absorption 



on uranium nuclei 
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FIG. 2: Energy spectra of a particles for the given angles and the given targets. 



III. DATA ANALYSIS 



Cross sections were analyzed in terms of a simple model assuming a moving source pre- 
scription. For completeness the content of the model is briefly repeated here. Suppose 
an IMF is emitted statistically from a source. The intensity distribution, by assumption 
a Maxwell-Boltzmann distribution, is isotropic in the rest system of the source. In the 
laboratory system we then have 



dVtde 



Cv^exp 



— ( e — y2mev cos 9 + -mf ^ J /T 



(4) 
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FIG. 3: Same as Fig. [21 but for ^He emission. 



where C is a normalization constant, 6 the emission angle and e the energy of the fragment. 
m is the mass of the fragment, v denotes the velocity of the source and T its temperature. 
In the present model v = v{e) and T = T(e) and not constants as in the usual moving source 
model. It is a common belief that in the early stage of a reaction the excitation energy is 
shared by a small number of nucleons. Thus momentum and energy conservation require 
a large source velocity and a high temperature in this stage, which is represented by the 
high energy of the IMF. At a later stage a succession of nucleon-nucleon interactions have 
taken place and more nucleons are in the source. This results in a smaller source velocity 
but higher temperature. How can one extract these two quantities? Unfortunately it is 
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FIG. 4: Same as Fig. |2but for ^Li emission. 



impossible. One can extract only a function of both quantities. The logarithm of the cross 
section is 



In 



dQde 



In(Cv^) + 



+ -mf ^ — \j2mtv cos 0^ /T 
= a(e) cos 9 + 6(e). 



(5) 



In the last line we have used the abbreviations 



6(e) = ln(Cv/^) - (^e + ^mv^^ /T 



(6) 
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FIG. 5: Same as Fig. |2lbut for ^Li emission. 



and 



a e 



V2 



niev 



T 



(7) 



We have chosen a and h in such a way to be consistent with earher nomenclature 
Linear fits to the logarithm of the double differential cross section versus the cosine of the 
emission angle are excellent. Both fit parameters a and h contain the source velocity and 
the temperature. Since h also contains the normalization constant it is impossible to deduce 
the numbers of interest. 

An emitted IMF is accelerated in the Coulomb field. In order to compare the energies 
before acceleration we study aj as function of e — Vc- This is done in Fig. El for the 
two targets aluminum and gold and for IMFs' a's, ^ife's, ^Lz's and ^Lz's. These are cases 
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FIG. 6: Same as Fig. I^lbut for ^Li emissfon. 



with reasonable or even good statistics. In the case of gold two components are visible: that 
for the higher energies is a smooth curve while for smaller energies a reflects the barrier 
penetration. It is interesting to note that the higher energy component follows almost a 
straight line with uniform slope. In order to show this effects we have fitted a straight line 
to the case of a-particle emission with the aluminum target. This line without any shift also 
passes through the bulk of points in case of other IMF types although not perfectly fitting 
the data. This may be an evidence that equilibration proceeds in an almost unique way 
independent of the target size. However, in order to extract angle integrated cross sections 
the fitted values with error bars were applied. 

In order to test the assumption that the low energy part is dominated by barrier pene- 
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FIG. 7: Same as Fig. [21 but for ^Be emission. 



tration we have fitted one single source with a constant temperature and a constant source 
velocity to the cross sections in the range from 15-30 MeV. The usual practice, as we have 
also applied above, namely, to correct the energy by subtracting the Coulomb barrier energy, 
is not applicable since it leads to negative energies. Consequently, Eq. (jH) can not reproduce 
the data. In this case one has to take the barrier penetration explicitly into account. We 
multiply the r.h.s. of Eq. (j3]) by the penetration probability 2^ 2^ 



Pie) 



huj 

27re 



In < 1 + exp 



27r (e - Vc 



c) 



(8) 



where uj is the frequency associated with a mean potential to be tunnelled through. Whereas 



fitting to an excitation function a+'^^^U fission 2&\ leads to /k^; ~ 4 MeV [2^, the present 
result is 6.0(3) MeV. The rather small Coulomb barrier of 17.67(23) MeV corresponds to 
a large radius of the emitting system. This might be an indication that the highly excited 
nucleus has expanded. For the source velocity the fit results in 0.0025(4)c, while one would 
expect 0.0033c from momentum conservation. This is an indication of fast particle emission 
in the equilibration process. The result of this exercise is shown in Fig. ^2 Finally we report 
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FIG. 8: Same as Fig. [21 but for ^Be emission. 



a fit value for the temperature of 3.12(17) MeV. There might be a correlation between the 
fit parameters. We have therefore performed moving source fits with barrier penetration 
to the angle integrated spectra, which will be discussed below. Unfortunately the different 
components are not clearly distinguishable as they are for the case of a-particle emission 
from gold, thus resulting in fits which are not so good. But again we find rather small values 
for the barrier. It will be interesting to study further data around the barrier and to see 
whether the barrier is reduced in comparison to a nucleus in its ground state. 

We use the slope and intercept parameters in Eq. (0) to get angle integrated cross 
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FIG. 9: Same as Fig. [21 but for ^''B emission. 



sections. The angle integrated cross section is 

^ = ^ {exp [6(e) + a(e)] - exp [6(e) - a(e)]} . (9) 

The resulting differential cross sections for the three targets are shown in Figs IT7HTII1 For 
the two lighter targets the energy distributions show an almost exponential slope without 
structure. In case of the gold target this structure is modified due to Coulomb effects. The 
distributions are discussed below. 



IV. DATA COMPARISON 



Although there are no data on IMF cross sections with exactly the same beam energy 
and for the same targets as employed in this study, there are data for energies or targets 
close by. We will compare the present data with those. First we will compare differential 
cross sections for the reaction p +^'' Al ^ {A = 7) + X . For that purpose the present cross 
sections for ''Li and ^Be emission were added. In Fig. these cross sections are compared 
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FIG. 10: The slope parameter as function of the Coulomb corrected energy. 



with those from Kwiatkowski et al. j22l| taken at a beam energy of 180 MeV. They measured 
fragments with masses A > Q and energies down to e/A > 0.05 MeV/u. The data have only 
a moderate overlap with the present data. There seems to be consistency between both data 
sets with respect to the absolute height as well as the shape of the spectra. 

As already stated in the introduction there are two studies of IMF emission from Ag. 
One was performed at a proton beam energy of 161 MeV Also these data cover 

smaller fragment energies than the present due to a gaseous AE detector. They observed 
fragments with charge number Z ranging from 3 to 12. The total cross sections from this 
measurement are shown in Fig. IT^a s function of the fragment charge number together with 
those from Green and Korteling [18| and the present results. For the case of the aluminum 
targets a large fraction of the cross section is missing due to the thickness of the first AE 
detector used here (see the latter comparison with model calculations and Figs. fT7lfT9|l . 
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FIG. 11: Fit with one source to the low energy data. The laboratory emission angle is given next 
to the appropriate data set or fit curve, respectively. The data and the fits were multiplied with 
the indicated factors. 



This makes the discrepancy in the yields esp. for Z = 5 and 6. The yield in the case of 
the cobalt target agrees best in that of the silver target. In Fig. IT^ angular distributions 
for Li and Be fragments integrated over the acceptance range are compared with those of 
Ref. Iisl. Again the agreement is reasonable with respect to the different energy ranges in 
the different experiments. Summarizing this comparison one can state that there is a fair 
agreement between the different measurements for Z = 3 and 4. It may be more instructive 
to continue the comparison on the level of spectra. 

This is done in Fig. ^Jfor the case of isotopic ^Li emission from Ag (210 MeV, Ref. [isj ]) 
and ^^Co. There is excellent agreement between the two data sets with respect to shape 
and absolute height. Finally we compare the slope of the present data in case of a-particle 



emission for the aluminum and gold target with those of ref. 14|. The later have been 
multiplied with an overall normalization factor of four. It becomes clear that the angular 
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FIG. 12: Double differential cross sections for A = 7 emission following the bombardment of 'Al 
with protons. The data from Ref. j^] are shown as histograms and are for 180 MeV, the present 
data are shown as dots with error bars. The emission angles are given in the legend where the first 
angle is for the present data, the second for those from Ref. [2^ . 



dependencies agree with each other in the overlap region. 

From these comparison it becomes evident that one can state that the present data are 
correct with respect to spectral shapes, angular distributions and absolute magnitude of the 
cross sections. 
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FIG. 13: Cross sections for total IMF production. The present data are shown by fuh symbols, 



those from Ref. by open dots and those by Ref. 18] by open squares. 
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FIG. 14: Angular distributions of Li and Be fragments for the indicated targets (this work) and 
for a silver target at 161 MeV proton energy 



V. COMPARISON WITH MODEL CALCULATIONS 

In this section model calculations are compared to the deduced data. A variety of models 
for fast particle emission in nuclear reactions are discussed in Ref. js^. A model especially 
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FIG. 15: Comparison of double differential cross sections from the p Co — >^ Li + X reaction 
(this work, dots with error bars) and p + Ag -^^ Li + X reaction (Eef . histograms). The first 
angle given is for the present data while the second is for those from 



suited for higher projectile energies is the intranuclear cascade model (INC). Although it 
cannot account for IMFs during the equilibration process it predicts the final excitation 
energy of an equilibrated system. This system then undergoes de-excitation by evaporating 
lower energy particles. 

The INC model was first proposed by Serber in 1947 [28^. The successful realization of 
this model by means of Monte Carlo simulations was published by Goldberger who did the 
first ca^ 
et al. 
INCL4 



culations by hand in 1948 |29[. Computer simulations were first done by Metropolis 



30 . In the present work we have applied the model in the standard Liege version 



The Intranuclear Cascade (INC) Model simulates - by the Monte Carlo method - se- 
quences of the nucleon-nucleon collisions proceeding inside the nucleus. This is equivalent 
to solving the Boltzmann transport equation for the time dependent distribution of the 
nucleons in the nucleus, treating explicitly collisions between the nucleons. As mentioned 
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FIG. 16: Similar as Fig. El but for the indicated reactions and data from Ref. jlj], normalized by 
a factor of four. 

above, such a picture of the reaction is justified in the case when the energy of the projectile 
is high enough. The INC is stopped when signatures are fulfilled, which indicate equilibra- 
tion of the decaying nucleus. In the INCL4.2 code the equilibration time r is determined by 
reaching a constant emission rate of cascade particles during the INC process. Typically r 
is of the order of lO^^^s or 30 fm/c. The longer this somewhat "artificially" chosen time the 
smaller E* is left for the evaporation process. Here we have chosen 

0.16 

/ \ 



208 



(10) 

but corresponds to the 



with To = 44.1 fm/c. This is smaller than the value used in 
earlier used value j33[. 

The description of each cascade involves three different stages: (i) initialization of the 
properties of the spatial and momentum distribution of nucleons in the nucleus, (ii) propa- 
gation of nucleons inside the nucleus, and (iii) collisions of the nucleons. 

The simplicity of the model and speed of calculations makes the INCL model very at- 
tractive. Of course the model cannot efficiently describe evaporation of the particles from 
the compound nucleus formed in the first stage of the reaction for two reasons: (a) the 
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evaporation is very sensitive to the density of states of the nuclei participating in the reac- 
tion, whereas the single particle density of states implicitly present in the INCL calculations 
is not exact enough, and (b) the calculations of the cascade over such long times as those 
characteristic for the compound nucleus emissions are not stable numerically and very ineffi- 
cient. The other very serious drawback of the INCL model is absence of correlations between 
nucleons, which could lead to emission of complex fragments. This is because the INCL is 
a single particle model with the mean field treated in an oversimplified manner. The mean 
field used in the INCL makes the assumption of being constant throughout the volume of 
the nucleus or modified at the surface of the nucleus, but it is always a static field. 

In practice the model of the intranuclear cascade (or equivalent) is applied to describe 
the first stage of the nuclear collision and the calculations are stopped once it can be as- 
sumed that equilibrium has been achieved. In the present study the INCL4.2 computer code 
was applied for this purpose. Discussion concerning the criteria for the terminating of the 
intranuclear cascade are presented by J. Cugnon et al. in Ref. [s^. 

After equilibration is reached we apply an evaporation model. It is the generalized evap- 



oration mode 
approach jsS, 



_(GEM) of Furihata |34|, which is based on the classical Weisskopf - Ewing 



36|. According to this approach, the probability of evaporation of the particle 



j from a parent compound nucleus i with a total kinetic energy in the center-of-mass system 
between e and e+de is defined as: 

Pj{e)de = gj(Tinv{e) ^'^^^ ,2\ — (11) 

where E is the excitation energy of the parent nucleus z, d denotes a daughter nucleus 
produced after the emission of ejectile j, and pi, pd are the level densities for the parent and 
daughter nucleus respectively. Q denotes the Q - value of the reaction. The statistical and 
normalization factor Qj is defined as gj = {2Sj + l)mj/7r^/i^, where Sj and mj are the spin 
and the mass of the emitted particle j respectively. The cross section ainv for the inverse 
reaction is evaluated from 

a^nvie) = agP{e)) (12) 

where cxg is the geometrical cross section. GEM considers fragments heavier than helium 
nuclei. There are 66 ejectiles (see Table H}. For the barrier penetration probability P we 
have used the form of Ref. |22i]. The parameters for light particles (n, p, d, t, ^He and ^He) 
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TABLE I: The ejectiles taken into consideration in the GEM calculations. 



Zj Ejectiles 
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8Li 


9Li 
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24Ne 
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21 Na 


22Na 


23Na 


24Na 


25Na 






12 


22Mg 


23 Mg 


24Mg 


25Mg 


26Mg 


27Mg 28Mg 



are taken from Ref. |^ whereas those for IMFs were adopted from the work of Matsuse et 
ah |38|. 

The total decay width Fj is calculated by integrating Eq. (fTTj) using Eq. (fT^ and is 

expressed as: 

r, = ^.f eP{e)p,{E -Q- e)de. (13) 

where V is the Coulomb barrier. For the level density we have applied the Fermi gas 
model expression 



2^a{E-S) 

p(E) = — — - ioTE>E^ (14) 

12 J. 

where a = A^/S (MeV~^) is the level density parameter, 5 is the pairing energy of the 
residual, T is again the nuclear temperature given by 1/T = \faJU~x — l.b/Ux where Ux is 
defined as Ux = 2.5 + IbQ/Ad. The excitation energy E^;, for which the formula for level 
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FIG. 17: (Color online) Energy differential cross sections for indicated IMFs for aluminum. The 
data are shown by the different symbols indicated in the figure. The histograms are calculations 
described in the text. 



density changes its form, is evaluated as Ex=Ux+S. To obtain a smooth continuity between 
the two formulae, the Eq parameter is determined as follows: 



Eo = E^-T -{InT - 0.25 In a - 1.25 In + 2 Vaf/^). (16) 

The contribution of the emission of IMFs in a long living excited state is taken into account 
together with those which decay to the ground state. The condition for the lifetime of 
excited nuclei considered in GEM is as follows: Ti/2/ln2 > h/V*. The value of F* is defined 
as the emission width of the decaying ejectile and is calculated in the same way as for the 
ground state, i.e. by Eq. (fT!?|) . The total emission width of an ejectile is summed over its 
ground state and all its excited states. All input parameters are the standard parameters of 
the models. We have not adjusted parameters to fit the present data. 

In Figs ITTHT^ we compare the angle integrated cross sections with the results of the 
calculations sketched above. A general trend is visible. For high IMF energies the calculation 
underestimates the experiment. The non-equilibrium fraction of the cross section is quite 
large in agreement with other experiments Q, The heavier the IMFs are, the more the 
agreement between the calculations and data deteriorates, even in the evaporative region. 
Heavy IMFs are strongly underestimated. In the case of gold there is emission of IMFs 
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from lighter composite systems, most probably fission fragments. This is also visible in the 
calculation but to a lesser extent than is observed in the data. The sequence within the 
isotopes is always obeyed by the calculations. However, in the case of '^Be emission from 
a gold-like composite the calculation predicts an almost negligible cross section, while the 
experiment is orders of magnitude higher. This is especially true for low energies and may 
point again to emission from lighter sources than the target-like system. 

VI. DISCUSSION 

We have measured IMF (He-C) emission for ^^Al, ^^Co and ^^^Au targets at a proton 
beam energy of 200 MeV, which is near the maximal abundance in the proton distribution 
in cosmic rays. The fragments were isotopically resolved. Spectra were taken at laboratory 
angles from 20°to 120°. Analysis in terms of a model of a moving source with continuous 
temperature and source velocity shows a linear relationship between these two quantities 
as a function of particle energy. The data in the case of gold show a strong influence of 
the Coulomb barrier. In the cases of the two lighter targets this feature was suppressed by 
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FIG. 19: (Color online) Same as Fig. El but for gold. 

the thickness of the first AE counter. Emission of fragments with a significantly smaller 
Coulomb barrier than for a target-like system is observed. The assumption that we observe 
emission from excited fission fragments was studied in evaporation calculations. Indeed, the 
calculations also show such fragments (see a-particle emission in Fig. [TU|) . although with a 
cross section more than one order of magnitude smaller than the experiment. The data were 
compared to model calculations. The first stage was calculated with an intranuclear cascade. 
In this cascade the emission of pions and nucleons can take place. After equilibrium is 
reached the energy and momentum distribution of the excited composite is transferred to an 
evaporation model, which, in addition to nucleons, allows for IMF emission. The frequency 
of isotopes being emitted for a specific element is followed by the calculation. The high 
energy tails visible in the experiments are not reproduced by the calculations. The emission 
of the heavier isotopes like ^Li, ^°Be, boron and carbon but also ^Be is underestimated. 
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The reduction of the Coulomb barrier observed may be due to dilution of the composite 
system. But an effect originating from the fission fragments can not be excluded. For the 
other two target nuclei we could not measure below the Coulomb barrier. Such data are 
highly desired to answer this question. Calculations, treating also IMF emission in the first 
fast stage are not yet available for the present data but are also desired. 

Let us now come back to the problem of Galactic cosmic rays (GCR) as discussed in the 
introduction. The production of BeB by Galactic cosmic-ray spallation of interstellar CNO 
nuclei was the standard model for BeB nucleosynthesis for almost two decades after first 
being proposed However, this simple model was challenged by the observations of 

BeB abundances in Population II stars, and particularly the BeB trends versus metallicity. 
Measurements showed that both Be and B vary roughly linearly with Fe, a "primary" scal- 
ing. In contrast, standard GCR nucleosynthesis predicts that BeB should be "secondary" 
versus spallation targets CNO, giving Be oc 02 Q. If O and Fe are coproduced (i.e., 
if the ratio 0/Fe is constant), then the data clearly contradict the canonical theory, i.e., 
BeB production via standard GCR's [J]. In order to accurately calculate the effects of the 
propagation of cosmic ray nuclei in the galaxy one needs to incorporate at least several 
hundred secondary cross sections into the propagation calculation. For charges with Z < 28 
his involves the fragmentation from ^ 55 nuclei with mass numbers A between 6 and 60 



42| |. The present data should help to improve our understanding of the systematics of the 



cross sections clS 8b function of Z, A, and A/Z. 

In GCR's one observes only stable isotopes since short lived isotopes decay. Thus only 
^He is observed because all tritium decays into it. The only difference might be ^^Be which 
has a half life of 1.6 x 10^ years. We therefore compare the ratio between ^^Be and ^Be for 
the different targets. The yield of the short lived ^Li is added to the latter. We report the 
ratio of the cross sections integrated over an energy range from 30 to 50 MeV in Table ITTl 
Within this range for all three particle types data exist. The ratio in cases of the two lighter 
targets is within error bars identical. For gold the primordial abundance of ^^Be relative to 
the A = 9 fragments is much larger than in the case of the two lighter targets. These ratios 
should be essential to study the age of GCR's. 
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TABLE II: The ratio R = a(}^Be)/[a{^ Be) + a{^Li)] for the different targets in the energy range 
30 to 50 MeV. 



target 


R 




0.654 ± 0.073 




0.577 ±0.133 




0.918 ±0.109 
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